Time periodic solutions of compressible fluid models of 

Korteweg type 

Zhengzheng Chen 

School of Mathematics and Statistics 
"^" I Wuhan University, Wuhan 430072, China 

Q . Qinghua Xiao 

CN I School of Mathematics and Statistics 

Wuhan University, Wuhan 430072, China 



S 



Huijiang Zhao* 

School of Mathematics and Statistics 



Abstract 



0^ 

CN ; Wuhan University, Wuhan 430072, China 

< 

C^ ' This paper is concerned with the existence, uniqueness and time- asymptotic stabihty of time 

periodic solutions to the compressible Navier-Stokes-Korteweg system effected by a time periodic 
external force in R" . Our analysis is based on a combination of the energy method and the time 
decay estimates of solutions to the linearized system. 
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O 
CN : 1 Introduction 

The compressible Navier-Stokes-Korteweg system for the density p > and velocity u = {ui,U2, • • • ,Un) 
k> ■ G M" is written as : 

_S_. (pt + V-ipu) = 0, 

\ {pu)t + V ■{pu<^u) + VP{p) - pAu - {u + p)V{V • u) = KpVAp + pf{t, x). 

Here, {t, x) G M+ x M", P = P{p) is the pressure, p, v are the viscosity coefficients, k is the capillary 
coefficient, and f(t,x) = {fi, f2, f3){t,x) is a given external force. System (1.1) can be used to 
describe the motion of the compressible isothermal fluids with capillarity effect of materials, see the 
pioneering work by Dunn and Serrin [1] , and also [21 [3l H] . 

In this paper, we consider the problem (1.1) for {p,u) around a constant state (poo)O) for n > 5, 
where poo is a positive constant. Throughout this paper, we make the following basic assumptions: 

(HI): p, V and k are positive constants and satisfying v -\ — p > 0. 
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(H2): P{p) is smooth in a neighborhood of poo satisfying P'{poo) > 0. 

(H3): / is time periodic with period T > 0. 

The main purpose of this paper is to show that the problem (1.1) admits a time periodic solution 
around the constant state {poQ,0) which has the same period as /. By combining the energy method 
and the optimal decay estimates of solutions to the linearized system, we prove the existence of a 
time periodic solution in some suitable function space. Notice that some similar results have been 
obtained for the compressible Navier-Stokes equations and Boltzmann equation, cf. [51 [9} lim fTT]. 

Precisely, Let A^ > n + 2 be a positive integer, define the solution space by 



Xm(0,T) 



{p,u){t,x) 



u{t, x) e C(0, T; /7^-i(M")) n C\0, T- H^-^{W^)), 

Vp{t,x) e L2(0,r;F^+i(M")), 

Vuit,x) e L^iO,T;H''iR^)),\\\ip,u)\\\ < M, 

for some positive constant Af and with the norm 



(1.2) 



lll(P,^)lll^ 



sup { 

0<t<T 



N + \\U 



-i}+[ (l|Vp(t)||^+i 
Jo 



+ Wu 



N, 



dt. 



(1.3) 



Then the existence of the time periodic solution can be stated as follows. 



Theorem 1.1. Let n > 5,N > n + 2. Assume the assumptions (H1)-(H3) hold, and f{t,x) £ 
C{0,T;H^-'^{W'-) n Li(R")). Then there exists a small constant 6o > and a constant Mq > 
which are dependent on p^, such that if 



sup ||/(i)||/fiv-inLi < ^0, 
o<t<r 

then the problem (1.1) admits a time periodic solution {f/"^^,u^'^^) with period T, satisfying 



(1.4) 



ipP' 



per\ 



eXM,io,T) 



Furthermore the periodic solution is unique in the following sense: if there is another time peri- 
odic solution {p^T^^ ,u^^^) satisfying (1.1) with the same f, and {ff^^ — poo,u^i^) G Xmo(0,T), then 

To study the stability of the time periodic solution (ff^^,u^^^) obtained in Theorem 1.1, we 
consider the problem (1.1) with the following initial date 



{p,u){t,x)\t=o = {po,uo){x) -^ (/9oo,0), as \x\ 



oo. 



(1.5) 



Here poix) and uo{x) is a small perturbation of the time periodic solution [pP^"^ ^u^^"^). And we have 
the following stability result. 

Theorem 1.2. Under the assumptions of Theorem 1.1, let {pP^'^,u^^^) be the time periodic solution 
thus obtained. If the initial date {po,uo) be such that \\{po — pP'^'^{0),uo — uP'^^{0)\\]\[-i is sufficiently 
small, then the Cauchy problem (1.1), HI. 5) has a unique classical solution {p,u) globally in time, 
which satisfies 

p_ pPer ^ C(0, oo; F^-i(E")) n C\0, oo; H^-^W)), 

u-uP^^ G C{0,oo;H^-^{W))nC\0,oo;H^-^{W)). 
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Moreover, there exists a constant Cq > such that 

Up - f^nml-, + \\{u- uPnm%-2 + J^ (I|V(p - ffn{r)\\%.^ + ||V(n - un{r)\\%.2) dr 

< Co (iipo - pp''m%-i + ho - uP'^m\%^2) , 

(1.7) 
for any t > and 

||(p-/>P^^u-nP^'■)IU- ^0 ast^oo. (1.8) 

Now we outline the main ingredients used in proving of our main results. For the proof of 
Theorem 1.1, thanks to the time decay estimates of solutions to the linear system ()2.7p (see Lemma 
2.1 below), we can show the integral in (|4.5|) is convergent. Based on this and the elaborate energy 
estimates given in Section 3, we prove the existence of time periodic solution by the contraction 
mapping principle. Here, similar to the case of compressible Navier-Stokes equations, Theorem 1.1 
is obtained only in the case n > 5 because of the convergence of the integral in ()4.5p . Thus, how to 
deal with the case n < 5, especially, the physical case n = 3, is still an open problem. Theorem 1.2 
is established by the energy method. The key ingredient in the proof of Theorem 1.2, among other 
things, is to get the a priori estimates, which can be done similarly to the estimates in Section 3. 

There have been a lot of studies on the mathematical theory of the compressible Navier-Stokes- 
Korteweg system. For example, Hattori and Li |12i [T3] proved the local existence and the global 
existence of smooth solutions in Sobolev space. Danchin and Desjardins [^ studied the existence 
of suitably smooth solutions in critical Besov space. Bresch, Desjardins and Lin [5] considered the 
global existence of weak solution, then Haspot improved their results in [B]. The local existence of 
strong solutions was proven in [T3] . Recently, Wang and Tan [15] established the optimal decay rates 
of global smooth solutions without external force. Li [16] discussed the global existence and optimal 
L^-decay rate of smooth solutions with potential external force. 

The rest of the paper is organized as follows. In Section 2, we will reformulate the problem and 
give some preliminaries for later use. In Section 3, we give the energy estimates on the linearized 
system (j2.4p . The proof of Theorem 1.1 is given in Section 4. In the last section, we will study the 
stability of the time periodic solution. 

Notations: Throughout this paper, for simplicity, we will omit the variables t, x of functions 
if it does not cauchy any confusion. C denotes a generic positive constant which may vary in 
different estimates. (•,•) is the inner product in Li^iW^). The norm in the usual Sobolev Space 
H^{W^) are denoted by || • \\s for s > 0. When s=0, we will simply use || • ||. Moreover, we denote 

n 

II • Wh- + II • llii by II • IIh^oli- If 5 = {91,92,- ■■ ,9n), then \\g\\ = ^(||fifcf )2. V = {di,d2,--- ,dn) 

k=l 

with di = dxi,i = 1, 2, • • • , n and for any integer / > 0, \7 g denotes all x derivatives of order / of the 
function g. Finally, for multi-index a = (ai, a2,- ■ ■ , ctn), it is standard that 



i=l 

2 Reformulated system and preliminaries 

We reformulate the system (1.1) in this section. Firstly, set 



Oi 



f^Tl T t Poo , II , V^ll . 7 . Poo 

7=V-P(Poo), K= K, fi = , U= , Al = , A2 = : 

7 Poo Poo Poo 7 
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and define the new variables 

CT = p-poo, V = X2U, 

then the system (1.1) is reformulated as 

( at + jV -v = Gi{a,v), 

\vt- fi'Av - u'\7{\7 ■ v) + -fVa - k'VAct = ^2(0-, v) + A2/, 

where 

Gi{a,v) = -XiV-{av), 



G-2{a,v) = -^ -{p/^v + vV {V ■ v)) -\i{v-V)v-\2 



(2.1) 



V(j. 



(2.2) 



Poo{(y + Poo) i Cr + Poo Poo 

Notice that Gi and G2 have the following properties: 

Gi((T, v) ~ Vo" • V + aV ■ V, 

G2{a, v) ~ (jAi; + o-V(V ■ v) + {v ■ V)v + aVa. 

Here ~ means that two side are of same order. 
Set U = {a,v),G= (Gi,G2), F = (0,A2/) and 

(0 "-^div \ 

7V - k'V A - ^' A - v'Vdiv J 

then the system (2.1) takes the form 

Ut + MJ = G{U) + F. (2.3) 

We first consider the linearized system of (|2.ip : 

\at + ^V-v = Gi{U), 

{ ~ (2-4) 

y vt - p'Av - u'V{V ■ v) + 7Va - n'VAa = G2{U) + A2/, 

for any given functions U = {&,v) satisfying 

a E F^+2(M"), V G F^+i(M"). 

Notice that the system ()2.4p can be written as 

Ut + AU = G{U) + F. (2.5) 

By the Duhamel's principle, the solution to the system ()2.4p can be written in the mild form as 



U{t) = 8{t,s)U{s)+ §{t,T){G{U) + F){T)dT, t>s, (2.6) 

J s 

where §(t, s) is the corresponding linearized solution operator defined by 

S(t,s) = e(*-^)*, t>s. 
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Indeed, the corresponding homogeneous linear system to (|2.4p is 

o"t + 7V • V = 0, 

vt - fi'Av - u'Viy ■ v) + 7Vcr - k'VAct = 0, 

a\t=s = crs{x), v\t=s = Vsix). 



(2.7) 



By repeating the argument in the proof of Theorem 1.3 in y^, we can get the following result for 
the problem (j2.7|) . The details are omitted here. 



Lemma 2.1. Let I > be an integer. Assume that {cr,v) is the solution of the problem (2.1) with 
the initial date Us G -ff'^^ Pi L^ and Vg G H^ Cl L^, then 



<Cil + t)"i\\ia„v,)\\Li + \\ias,Vs)\\), 

IIV^+MOII < C(l + t)-t-^ {\\{cTs,Vs)\\l^ + IKV^+V,, V'^^,)! 

\\V'vit)\\ < C(l + t)-t-t (||(a„^.)||z,i + ||(V^'+V„ V'^t;,)!!) 
where k is an integer satisfying < k < I. 



3 Energy estimates 

In this section, we will perform some energy estimates on solutions {a, v) to problem (|2.4p . Through- 
out of this section, we assume that f{t,x) G H^~^{W^) n L^{W^) for all t > 0. For later use, we list 
some standard inequalities as follows, cf. [8]. 

Lemma 3.1. Let m be a positive integer and u G ff'2'^ ^W"), then 

\\u\\loo < CllV^+^nllllV^-^nll for n = 2m, 

\\u\\loo < CIlV^+^-ull IIV^ull for n = 2m + 1. 

Lemma 3.2. Let m be the integer defined in Lemma 3.1 and f,g,h G H'-^''^ (M") , then we have 



(^) 



in) 



f ■ g ■ hdx 



<e||V™-Vll^ + a||5l 



< e 



I r< ||v™~-^nlrll/)lr 



f ■ g ■ hdx 

for any e > 0. Here and hereafter, C^ denotes a positive constant depending only on e. 

We first give the energy estimate on the low order derivatives of {(J,v). 

Lemma 3.3. Let n > 5, N > n + 2, then there exists two suitably small constants do > and eo > 
such that for < e < eo, it holds 



d 



- {\\U{t)f + \\Va{t)f + do{v,Va){t)) + \\Vv{t)f + ||Va(t)||? 
where m is defined in Lemma 3.1 and C depends only on p^,^,v and k. 



(3.1) 



(3.3) 
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Proof. Multiplying (j2.4p i and ()2.4|) 9 by a and v, respectively, and integrating them over M", we 
have from integrating by parts that 

= {Gi{U),a) + {G2{U),v)+K'{VAa,v) + \2{f,v) (3-2) 

= /0+/l+/2 + /3- 

From (|2.2|) and Lemma 3.2, we have 

lo < e|| V^-Vlli + aC (llVaf ||i}||2 + ||af ||Vt}f ) 
< e||V™-V||i + C,C||f/||2||Vi/||2, 

and 

h < ellV^-^^lli + C7,C7||{7f||VC/||?. (3.4) 

For I2, integrating by parts and using l\2.^ ^ . ()2.2p and Lemma 3.2, we deduce that 

l2 = -K'{Aa,V-v) = -{Aa,at-Gi{U)) 

7 

= -f^\\Var--{Aa,G,m (3.5) 

2^ at 7 

— rll^'^^ll +^l|v '^11 + L-eO II vty II ||v ^yib- 

27 at 

For /3, Lemma 3.1 gives 

l3<e\\V'^-'v\\l + C,C\\f\\l. (3.6) 

Since n>5, A^>n + 2, we have m — 1 > 1. Substituting (|3.3|) - (|3.6|) into (|3.2|) yields 



^^{\\Uf + \\Va\\') + \\Vv\\' + \\V-vf 

< eC (II V™- V||i + llV^af) + 6C7||V2^||^_i + C,C {\\U\\l+i l|VC/||? + 



(3.7) 



Li ' 



provided that e is small enough, where G depends only on p^o, /^, z^ and k. 

Next, we estimate ||V(t|P. Taking the L^ inner product with Vo" on both side of ()2.4p 9 and then 
integrating by parts, we have 

7||V(Tf + K'||VVf 

= -{vt,Va) + fi'{Av, Va) + iy'{V{V ■ v), Va) + (^2(^7) + A2/, Va) (3.8) 

= h + h + h + h- 
Similar to (3.5), the term /4 can be controlled by 

h = -j^{vy(j)-{V -v^at) 

^'Jt^""' ^""^ - (V • «, -7V • « + Gi(c7)) (3.9) 

< -j{v,Va) + 27IIV • vf + C||V™-i^||^||V[7f . 
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Integrating by parts and using the Cauchy-Schwartz inequality, it is easy to get 

h + h< jll VVf + C{\\Vvf + ||V • vf). (3.10) 

Finally, (|2.2p and the Cauchy-Schwartz inequality imply that 

h < |||Vaf + C {\\V"'"^U\\l\\VU\\j + ll/f ) . (3.11) 

Combining ()3.8p - (|3.1ip . we obtain 

^{v,Va) + \\Vaf + \\V'af 
at 



< C{\\Vvf + ||V • wf) + C \\V"'-Wg\\VU\\l + 



(3.12) 

where the constant C depends only on Poo, fJ-, ^ and k. Multiplying ()3.12p with a small constant do > 
and then adding the resultant equation to (|3.7p . one can get (|3.ip immediately by the smallness of 
do and e. This completes the proof of Lemma 3.3. 

Next, we derive the energy estimate on the high order derivatives of {(J^v). We establish the 
following lemma. 

Lemma 3.4. Let n>'a,N>n + 2, then there exists two suitably small constants di > and ei > 
such that for < e < ei, it holds 

Y, I llVa(t)||^ + ||Vt;(t)||^_i + di f;(a>,a:Va)(t) I +||VMi)ll^ + l|VMi)ll^-i , , 

^^ \ \a\=l J (3-13) 

2 I r^ r< fllV7T~TU\ll2 llV7T~Tr4-\ll2 , ll J^/'+M|2 



< eC\\Va{t)r + aC (||VC/(t)||^_2l|Vf/(t)||^ + ||/(t)||^-ij , 
where C is depending only on poo,fM,i' and k. 

Proof. For each multi-index a with 1 < |a| < N, applying d^ to ()2.4p i and ()2.4p 9 and then taking 
the L^ inner product with d^a and d^v on the two resultant equations respectively, we have from 
integrating by parts that 

1 d 



\d>f + ||9>f ) + p-'\\d^Vvf + i^'Wd^V ■ v\ 



111// cjii -i- \\i / 7 711 I -t- / /, 1 1 ^ / \/ 7 ; 1 1 -t- ;y 1 1 / / \/ • / ; i 

2dt 

{d^G,{U),dSa) + {dSG2{U),dSv) + K'{dSVAa,d^v) + X2{dSf,d^v) (3-14) 



= /8+/9 + /l0+/ll- 

Now, we estimate Is-Iii term by term. For /§, we deduce from (j2.2p and the Cauchy-Schwartz 
inequality that 

(3.15) 
< eWd^af + aC {\\d-{Va ■ v)f + \\d-{aV ■ v)p) . 

By Leibniz's formula and Minkowski's inequality, we get 

0<m = \a\-l 

+C Y. C^Wd^.Va-d^^vf 

0<|/3|<|a|-2,|a-/3|<f (3.16) 



+C Y. C^\\d^.Va-dr 



'vP 



0<|/3|<|a|-2, |a-/3|>^ 



2 



J0 + J1 + J2 + h- 
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Here C? denotes the binomial coefficients corresponding to multi-indices. For Jq, lemma 3.1 gives 

Jo < C (||^||io.||5-Vcr||2 + ||Vc.||2^||9a5||2) 

< C {\\Vv\\l_,\\V'~a\\l_, + ||V2a||^_,||V^||^_,) , 

where, in the last inequality of ()3.17p . we have used the fact that m — 1 > 1 and m + \ < N — 4 due 
to A^ > n + 2 and n > 5. Similarly, it holds that 

Ji<C Yl ll^r^^lli->l|5fVa|l' < C\\VM%-5\N'^\\%-2- (3.18) 

0<|^| = |q|-1 

For the terms J2 and J3, notice that for any /3 < a with |a — /3| < y, 

N n N N 

\a-f3\+m + !<- + - + !<- + - = N, 

and for any j3 < a with |a — /3| > -^, 

N n 
|/3| + m + 2 = |a|-|a-/3|+m + 2<iV \ h2<iV + l. 

which implies |/3| + m + 2 < A^ since |/3| and m are positive integers. Hence, we deduce from Lemma 
3.1 that 

J2<C Yl \\dr^v\\l^\\dlVaf < C\\V'v\\%.2\N'a\\%_,, (3.19) 

0<|/3|<|a|-2, |a-^|<f 

and 

Js<C Y ||5fVa||ioo||ar^t)f < C\\VM%_,\\V'a\\%_,. (3.20) 

0<|/3|<|a|-2,|a-^|>f 

Putting (f3T7D - (f3:20D into (I3J6I1 . we arrive at 

||5^(Va • v)f < C (||V^||^_5||V2a||^_i + ||V2c7||^_3||VC7||^_i) . (3.21) 

Similarly, it holds 

miaV ■ v)f < C [\\Va\\%_,\\V^v\\%_, + ||V2[/||^_3||Vt/||^_i) . (3.22) 

Combining (|3J5]) . (|3:2T]) and ([3:22]) yields 

h < 4d>f + <^eC||Vt/||^_2||Vt/||^. (3.23) 

For the term /g, let oq < a with |ao| = 1) then 

I9 = -(9°-°oG2,9°+°°i;) < e||5°+°0i;f + Ce||9°-°°G2f . (3.24) 

Similar to the estimate of ()3.2ip . we have 

||9°-°°G2f < C||[/|||^_i||VC/||^. (3.25) 
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Thus, it follows from (l3:2i]l and ([3:25|) that 

h<e\\d'^+^^^vf + C\\U\\l_,\\VU\\l. (3.26) 

Notice that ^^?n\ and (f3:22]l imply 

Wd'^G.f < C iWd^iVa ■ v)f + ||a:(aV • v)f) < C7||V^||^.2l|V^||^. (3.27) 

Therefore, we derive from ()2.4p i . p.27p and the Cauchy-Schwartz inequality that 

7 

= -^{d^Va,d^Vat) - -{d^Aa,d^Gi{U)) (3.28) 

7^ 7 

- "1^1"^"^^"' + 4d^^af + C,C7||VC7||^_2l|Vf/||^. 
Moreover, it holds that 

111 = -A2(9r°°^,5r"''/) < e||9r"°Hl' + aC\\f\\%_,. (3.29) 

where oq is defined in (l3:2iD . Combining UiSAM . dS^SD, dS^S]), (l3:28|) and IK29\\ . if e is small enough, 
we have 

I {\\d>\\l + ||5>in + ||5,"Vt;f + \\d^V ■ vf 

< eC\\dyr + eCWd^Aa^ + ^C (||V?7||^_2l|Vf/||^ + ||/||^„i) , 

where C depends only on Poo, fJ-, t^ and k. 

Now we turn to estimate ||(?"A(t|P for 1 < |a| < A^. As we did for the first order derivative 
estimate, applying d% to ()2.4p 9 and then taking the I? inner product with d^Va on the resultant 
equation, we get from integrating by parts that 

K'Wd^Aaf + ^Wd^Vaf 

= -{d^vt,d;^Va)+p'{d^Av,d;^Va) + u'{d^V{V ■v),d^Va) 
+ {dSG2{U),d^Va) + X2{d^f,d^Va) 

= /l2+/l3 + /l4 + /l5+/l6- 

The first term /12 is controlled by 

h2 = -j^{d^v,d^Va) + {dy,d^Vat) 

= ~{d^v, d^Va) - {d^V ■ V, d^{-^V ■ V + Gi(c7))) (3-32) 

< -^(5>, d'^Va) + 27||a-V • vf + C||V?7||^_2||V^||^. 

Here, in the last inequality of ()3.32p . we have used ()3.27p . By integrating by parts, the Cauchy- 
Schwartz inequality and ()3.25p . the other terms /13-/15 can be estimated as follows. 

/i3 + hA < yII^x Wf + G {Wd^Vvf + \\d^V ■ vf) , (3.33) 



(3.31) 
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/i5 < j||a:+"«VcTf + C||[/||^_i||VC/||^, (3.34) 

/i6< jll^r'^^Vaf + C||/||^_i. (3.35) 

where ao is given in (j3.24p . Combining (|3.3ip - (j3.35p . we obtain 



< c {Wd-Vvf + \\d-v ■ vf) + c (||f/||^_il|vc/||^ + \\f\\%_,) . 



(3.36) 



Multiplying ()3.36p with a suitably small constant di > and then adding the resultant equation to 
(I330]l gives 

^ {\\d>\\\ + ||5>f + di(5>, d'^Va)) + Wd'^VaWl + ||9,"Vt; 
<eC\\d'^af + CCA\\VU\\%_^\\VU^^'^ ' n^n^ 



dt ^" ^ "^ II x II -,..... ,, ,, .. ,,. ,, .. ,, {2,.2>1) 



provided that di and e are small enough, where C depends only on p^,^,v and k. Summing up a 
with 1 < |a| < A^ in (|3.37p . then p.l3p follows immediately by the smallness of e. This completes 
the proof of Lemma 3.4. 

As a consequence of Lemmas 3.3-3.4, we have the following Corollary. 

Corollary 3.1. Let n > 5, N > n + 2, then there exists two suitably small constants do > and 
di > such that 

I (lk(t)ll^+i + IKt)ll^ + 4(^,v<T)(t) + di f;(a>,a,"v^)(t)j +||Va(t)||^+i + l|vKOI^^ 

|2 IIV7r7/+Ml2 I II ^^+M|2 \ 

(3.38) 



< C (||C/(t)||^_il|VC/(t)||^ + \\f{t)uHN-.nL^j , 

where C depends only on p^,^,v and k. 

Proof. Notice that, from the fact that ?n, — 1 > 1 and m + 1 < A^ — 4, we have 

\W^'yfm-2 + \WML-i<c\\v^u\\l_,, 

and 

WTTW ^ 1 1 7T 1 1 2 

W-J \\m+l ^ W-J llAf-4- 

Adding ()3.37p to ()3.ip . we obtain ()3.38p immediately by the smallness of e. This completes the proof 
of Corollary 3.1. 

4 Existence of time periodic solution 

In this section, we will combine the linearized decay estimate Lemma 2.1 with the energy estimates 
Corollary 3.1 to show the existence of time periodic solution to (1.1). Now, we are ready to prove 
Theorem 1.1 as follows. 
Proof of Theorem 1.1. The proof is divided into two steps. 

Step 1. Suppose that there exists a time periodic solution U^^'^{t) := {a^'^'^ {x , t) , v^^'^ {x , t)) , t G M 
of the system ()2.ip with period T, and [/^^^(t) G Xmq{0,T) for some constant Mq > 0. Then it 
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solves (|2.3p with initial date Us = U^^'^{s) for any given time s E M. Choosing s = —kT for A; G N. 
Clearly, UP'^^{—kT) = U^'^^{0), thus (|'2.3|) can be written in the mild form as 

[/P^'^(t) = S(t,-A:T)C/P"'^(0)+ / S{t,T){G{UP''''){T) + F{T))dT. (4.1) 

J-kT 

Denote S{t,-kT)UP^''{0) := (af ^(t), wf''(t)). Applying Lemma 2.1 to S(t, -A:T)C/p^'^(0), we have 

\wrmN <{i+t+ kT)-^. (ii(ar,^r)iiLi + ikrii^ + ii^ni^-o ,, ^, 

— > as K — )• oo. 
and 

\\vrmN-i < (1 + 1 + fcT)-? (ii(ar,^'r)iiL^ + ii^rii^ + ibni^-i) ^^ 3^ 

— > as k ^ 00. 

Since L^nL^ is dense in L^, g2|) and gSD still hold for UP^''{0) = ((Jg'^ v^"'') e ii'^(M'^) xi/^-i(M'^). 
On the other hand, denote 

S{t,T){G{UPn{r) + F{t)) := {Si{t,T),S2{t,T)). 

By using Lemma 2.1 again, we get 

||5i(t,r)|U<(l + t-r)-?iCo, ||S2(t,T)||jv-i<(l + t-T)-5i^o, (4.4) 

where 

Ko = ||(Gi([/P-),G2([/P^n + A2/)(r)|ki 

+ ||Gi(C/P-)(r)|U + \\iG2iUPn + A2/)(t)||^_i. 

Then (|4.4p guarantees the convergence of the integral in ()4.ip since f > 1 when n > 5. Thus, letting 
/c — )• 00 in (|4.ip . we obtain 

UP^'-{t)= f S(t,r)(G(C/P^'^) + F)(r)(ir. (4.5) 

./ — 00 

For any ?7 = {a,v) G XAfo(0,T), define 

num= [ S(t,r)(G([/)+F)(r)dr. 



—00 



Then (03]) shows that t/P^'" is a fixed point of ^[[/]. 

Conversely, suppose that ^ has a unique fixed point, denoted by Ui{t) = {ai,vi){t). We show 
that Ui{t) is time periodic with period T. To this end, setting U2{t) = Ui{t + T). Since the period 
of / is T, the period of F is T too. Thus, we have 

U2{t) = Ui{t + T) = ^Ui]{t + T) 

S(t + T,r)(G([/i)(r)+F(r))dr 

00 

/ Sit + T,s + T){GiUi)is + T) + Fis + T))ds (4-6) 

t 

S{t,s)iG{U2){s) + F{s))ds 

-00 

M/[[/2](t) 
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where we have used 

§(t + r,s + r) = s{t,s). 

Then by uniqueness, U2 = Ui, which proves the periodicity of Ui{t). Since Ui{t) is differentiable 
with respect to t, it is the desired periodic solution of the system (2.1). 
Step 2. Now, it remains to show that if (H1)-(H3) hold, and 

sup ||/(i)||i/iv-inLi 

0<i<T 

is sufficiently small, then ^ has a unique fixed point in the space Xj\/p(0,T) for some appropriate 
constant Mq > 0. The proof is divided into two parts. 

(i) Assume that U = {d',v) in the system (|2.4p is time periodic with period T. Denote U = ^[U] 
with U = {cr,v). Then by the same argument as ()4.6p . one can show that U is also time periodic 
with period T. Notice that U satisfies the system (j2.4|) . Thus, for n > 5 and N > n + 2, Corollary 
3.1 holds. Integrating ()3.38p in t over [0,T] to get 

(l|VcT(t)||^+i + ||V^;(t)||^)dt 


<c r (iif/(t)ii^-iiivc/(t)ii^ + \\fm%-i + wfmh) dt 

Jo ^ ' (4.7) 

<C7sup \\Um%-i [ \\VU{t)\\Ut+ I \\f{t)\\lN-,^L^dt 

0<t<T Jo Jo 

<C7|||c7"(t)ir + cr sup ||/(t)||^^-inLi- 

0<t<T 

On the other hand, by Lemma 2.1, we have 

\W{t)y< (l + t-r)-?i^i(ir, \\v{t)y-i< {l + t-T)-^KidT, (4.8) 

J ~oo J —00 

(4.9) 



where 

Ki = \\iGiiU),G2iU) + X2f)ir)hi 



+ \\Gi{U){t)\\m + \\{G2{U) + A2/)(t)|U_i. 
From (122]), ([3:25]) and (f3:27D . we easily deduce that 

||(Gi(c7)(r)||^i<C||VC/(T)||||?7(r)||, 
\\{G,{U){T)y < C||VC7(r)|U_2||V?7(r)|U, 
\\G2{U) + X2f){r)h^ < C||VC7(T)||ir(r)|| + C||/(r)||^i, 

\\G2iu) + X2f)ir)\\N-i < c\\u{T)y-i\\vuiT)y + c||/(t)|U-i 

Combining (|4.8p - (|4.10p . we obtain 



k(t)|U < C / (1 + t - r)-t (\\U{T)y-i\\VU{T)\\N + \\f{T)\\HN^^nL^) dr 
J —00 

<CY^A,+C {l + t-T)-^\f{r)\\HN^.^L-dT 

,_n J —00 



j=0 

00 



<CY,^j+C sup WimHN-inL^, 



(4.10) 



(4.11) 



,=0 o<t<T 
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where 



(4.12) 



Aj=C [ ' {l + t-T)-l\\U{T)\\N-l\\VUiT)\\NdT 

Jt~{j+1)T 

<C[ (i + t_r)-tdr / ||^7(r)||i._i||VC/(r)||^dr 

\Jt-{j+l)T j \Jt-{j+l)T J 

<c(i+ir)-t sup \\u{T)y-i(f \\vu{T)\\%dTy 

0<T<T \Jo J 

<c(i+ir)-f|||[/|||2 

Since \> \ when n > 5, substituting (|4.12p into ()4.1ip gives 

lk(t)IU<C|||C7|||2 + C sup ||/(t)||Hiv-inLi- (4.13) 

o<t<r 

Similarly, it holds that 

\\v{t)\\N-x<C\\\lJ\\f^C sup ||/(t)||H--inL- (4.14) 

0<t<T 

Thus, we deduce from (g^]), (1^^131) and (liJi]) that 

|||^[f/]|||<Ci|||C/|||2 + C72 sup ||/(t)|biv-inLi, (4.15) 

o<i<r 

where C\ and C2 are some positive constants depending only on poo. A*, v-, 1^ and T. 

(ii) Let C/i = (cTi,!;!) and C/2 = {^21^2) be time periodic functions with period T in the space 
-^Mo(Oi^)) where Mq > will be determined below. Then similar to (i), we can get 

WMUi] - nU2]\\\ < C-i (|||?7i||| + |||?72|||) \\\Ui - C/2III, (4.16) 

where C3 is a positive constant depending only on Poo; M; t^i ^ and T. Choose Mq > and a sufficiently 
small constant 5 > such that 

CiM^ + C72(5 < Mo, and 2C^Mq < 1 (4.17) 

That is, 

l-Vl-4CiC26 ^ • [ l + Vl-4CiC2^ 1 \ 

2C^ <Mo<mm| — ,— ,l| (4.18) 

Notice that 

1 - yi - 4CiC2<5 ^ ^ ^^^ 
^c^^ ^0 - '^^o. 

Then there exists a constant (5o > depending only on poo-,^J'-,v-,i^ and T such that if < (5 < (^Oi 
the set of Mq that satisfying ()4.18p is not empty. For < (5 < (5o, when Mq satisfies ()4.18p . ^ is 
a contraction map in the complete space Xa/q(0,T), thus ^ has a unique fixed point in Xmq{^-,T). 
This completes the proof of Theorem 1.1. 
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5 Stability of time periodic solution 

This section is devoted to proving Theorem 1.2 on the stabihty of the obtained time periodic solution. 
We shah estabhsh the global existence of smooth solutions to the Cauchy problem (1.1), (1-5). 

First, let {fp"^"^ ,vP^'^) be the time periodic solution obtained in Theorem 1.1 and (p, ti) be the 
solution of the Cauchy problem (1.1), (1-5). Denote 

{a,v) = {p- P00A2U). 
Let {d-,v) = {a — a^^'^ ,v — t;^^^), then {a,v) satisfies 

\vt- p'Av - u'V{V • v) + -fVa - K'VAa = G2(a + aP^\ v + yP^^) - G2(aP^^ yP^'' 

with the initial date 

a\t=o = cTo(x) = poix) - /^^(O), v\t=o = yo{x) = HM^) " ^^^'"(0)). 

Define the solution space by X(0,oo), where for < ti < ^2 ^ cO; 

a{t, x) £ C{h,t2; H^-HW)) n ^^^1,^2; //^-=^(M")), 
v{t, x) £ C{tut2; H^-^W)) n C\h,t2; H^-^W)), 
Va{t,x) G L'^{ti,t2;H^-~^{R'')),Vv{t,x) G L'^{ti,t2; H^''^{R'')), ^ 



(5.1) 



(5.2) 



X{h,t2) 



{a,v){t,x) 



} (5.3) 



with the norm 



{a,v){t)f:= sup {\\a 

ti<t<t2 



N-1 + \\V 



-2}+ Tdiv^ 
Jti 



N- 



-1 + l|Vw 



N~2, 



dt. 



(5.4) 



Notice that (ctP^^^p^'^) G X(0,r). 

By using the dual argument and iteration technique as |12] . one can prove the following local 
existence of the Cauchy problem (j5.ip , (|5.2p . We omit the proof here for brevity. 



Lemma 5.1. (Local existence) Under the assumptions of Theorem 1.1, suppose that (o^OjiJo) ^ 
// (M") X ff (M") and infpo(2;) > 0. Then there exists a positive constant Tq depending 
only on \\{ao,vo) j(\ such that the Cauchy problem /i5.1]) . /i5.2\) admits a unique classical solution 
(a,u) G X{0,Tq) which satisfies 

\\{a,vmS^<C4{ao,vo)l 

where C4 is a positive constant independent of ||((To,vo) -[f. 

As usual, the global existence will be obtained by a combination of the local existence result 
Lemma 5.1 and the a priori estimate below. 

Lemma 5.2. (A priori estimate) Under the assumptions of Lemma 5.1, suppose that the Cauchy 
problem i5.1\) . i5.2\) has a unique classical solution {(7,v) G X{0,Ti) for some positive constant Ti. 
Then there exists two small constants 5 > and C5 > which are independent ofTi such that if 



sup 
0<t<Ti 



{a,v)it))^<6, 



(5.5) 
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it holds that 

Mm%-i + ll^(*)ll^-2 + [ (l|Va(T)||^_i + ||Vt;(r)||^_2) dr < C, {\\ao\\%-i + \\vo\\%-2) (5-6) 

Jo 

for all te [0,Ti]. 

Proof. Noticing that some smallness conditions can be imposed on (a^^^jV^'^^), without loss of 
generality, we may assume |||(<7^*^'',u^^^)||| < e with e > being sufficiently small. Then by the 
similar argument as in the proof of Lemmas 3.3-3.4, we can obtain 

^ {\\Uf + \\Vaf + d2{v, Va)) + \\Vvf + \\Va\\l ^^ ^^ 

< eC {\\VM%-7 + W'^'nl-e) ^ 

and 

d ( ^~^ 

^ II Va||^_2 + II Vt;||^_3 + ^3 ^ {d'^v, d^Va) J + \\V'a\\%_, + \\VM%-3 (5,3) 

<eC{\\Vaf + \\Vvf), 

where d2 > and da > are some suitably small constants, and C is a constant depending only on 
Poo,fJ;i^ and K. Adding (j5.8p to (j5.7|) . it holds 

d ( ^"^ 

^ ll^ll^-i + 11^11^-2 + ^2(t^, Va) + d3 E^^-^'^"^^)j (5.9) 

+ ||Va||2^_i + ||V7;||^_2 < 0, 

provided that e is sufficiently small. Integrating ()5.9p in t over (0, t), one can immediately get (|5.6p 
since 

Af-2 

ll^ll^-i + l|i^ll^-2 + '^2(^,V(j) +(i3 Y, (5x^,5"Va) ~ ||a||^_i + ||V?;||^_2. 

|a|=l 

by the smallness of ^2 and ^3. This completes the proof of Lemma 5.2. 

Proof of Theorem 1.2. By Lemmas 5.1-5.2 and the continuity argument, the Cauchy problem 
(j5.ip . (|5.2p admits a unique solution {(t,v) globally in time, which satisfies (|1.6p and (|1.7p . Then all 
the statements in Theorem 1.2 follow immediately. This completes the proof of Theorem 1.2. 
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